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Abstract — In this work, the performance analysis of a dual- 
hop relay transmission system composed of asymmetric radio- 
frequency (RF)/free-space optical (FSO) links with pointing er- 
rors is presented. More specifically, we build on the system model 
presented in [1] to derive new exact closed-form expressions for 
the cumulative distribution function, probability density function, 
moment generating function, and moments of the end-to-end 
signal-to-noise ratio in terms of the Meijer's G function. We 
then capitalize on these results to offer new exact closed-form 
expressions for the higher-order amount of fading, average 
error rate for binary and M-ary modulation schemes, and the 
ergodic capacity, all in terms of Meijer's G functions. Our new 
analytical results were also verified via computer-based Monte- 
Carlo simulation results. 

Index Terms — Asymmetric dual-hop relay system, pointing 
errors, mixed RF/FSO systems. 

I. Introduction 

IN recent times, free-space optical (FSO) or optical wireless 
communication systems have gained an increasing interest 
due to its various characteristics including higher bandwidth 
and higher capacity compared to the traditional radio fre- 
quency (RF) communication systems. In addition, FSO links 
are license-free and hence are cost-effective relative to the 
traditional RF links. These features of FSO communication 
systems potentially enable solving the issues that the RF 
communication systems face due to the expensive and scarce 
spectrum [2]— [6]. However, the atmospheric turbulence may 
lead to a significant degradation in the performance of the FSO 
communication systems [2]. Additionally, thermal expansion, 
dynamic wind loads, and weak earthquakes result in the 
building sway phenomenon that causes vibration of the trans- 
mitter beam leading to a misalignment between transmitter and 
receiver known as pointing error. These pointing errors may 
lead to significant performance degradation and are a serious 
issue in urban areas, where the FSO equipments are placed on 
high-rise buildings [7]-[9]. 

On the other hand, relaying technology has gained enormous 
attention for quite a while now since it not only provides wider 
and energy-efficient coverage but also increased capacity in 
the wireless communication systems. As such many efforts 
have been made to study the relay system performance under 
various fading conditions [10]— [13]. These independent studies 
consider symmetric channel conditions i.e. the links at the 
hops are similar in terms of the fading distributions though 
it is more practical to experience different/asymmetric link 
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conditions at different hops i.e. each link may differ in the 
channel conditions from the other link [1], [14]— [17]. This is 
due to the fact that the signals on each hop are transmitted 
either via different communications systems or the signals 
might have to commute through physically different paths. For 
instance, as proposed in [1], a relaying system based on both 
FSO as well as RF characteristics can be expected to be more 
adaptive and constitute an effective communication system in 
a real-life environment. 

The model utilized in our work is similar to the one 
presented in [1]. Although [1] lacks the motivation behind 
such a model, we understand and proceed with such a model 
based on the following explanation. Considering an uplink 
scenario, besides all the advantages of FSO over RF, that very 
much motivates this work is the concept of multiplexing i.e. 
we can multiplex users with RF only capability into a single 
FSO link. This comes with the reasoning that there exists 
a connectivity gap between the backbone network and the 
last-mile access network and hence this last mile connectivity 
can be delivered via high-speed FSO links [18]. For instance, 
in developing countries where there might not be much of 
a fiber optic structure and hence to increase its reach and 
bandwidth to the last mile, it will require huge amount of 
economic resources to dig up the current brown-field. It 
will be much better to simply install FSO transmitters and 
detectors on the high-rise buildings and cover the last mile 
by having the users with RF capability to communicate via 
their respective RF bands and let the rest be taken care of 
by the FSO links to get it through to the backbone as can 
be observed from Fig. 1. This multiplexing feature will avoid 
the bottleneck situation for the system capacity and in fact be 
a faster option relative to traditional RF-RF communications 
wherein multiple RF's being sent through a single FSO at 
once. Hence, there are two outstanding features, among others, 
of this system to make it very advantageous over the current 
traditional system. Firstly, maximum possible RF messages 
can be aggregated into a single FSO link thereby utilizing 
the system to the maximum possible capacity. Simultaneously, 
the system benefits from another feature of having the RF 
link always available irrespective of FSO transmission since 
the RF and the FSO operate on completely different sets 
of frequencies allowing for no interference between them 
at any instant. Therefore, the RF frequency bands can be 
utilized by other possible devices/users around in range to their 
benefit while the FSO link is yet under operation. Above all, 
having FSO will avoid any sort of interference(s) also due 
to its point-to-point transmission feature unlike RF where the 
transmission is a broadcast leading to possible interference(s). 
In Fig. 1, there exists no fiber optics structure between the 
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ability (OP). In this work, we build on the model presented in 
[1] to study the impact of pointing errors on the performance 
of asymmetric RF/FSO dal-hop transmission systems with 
fixed gain relays. For instance, we derive the CDF, probability 
density function (PDF), moment generating function (MGF), 
and moments of the end-to-end signal-to-noise ratio (SNR) of 
such systems. We then apply this statistical characterization 
of the SNR to derive closed-form expressions of the higher- 
order amount of fading (AF), average bit-error rate (BER) of 
binary modulation schemes, average symbol error rate (SER) 
of Af-ary amplitude modulation (M-AM), M-ary phase shift 
keying (M-PSK) and A/-ary quadrature amplitude modulation 
(M-QAM), and the ergodic capacity in terms of Meijer's G 
functions. 



II. Channel and System Models 



Fig. 1, System model block diagram of an asymmetric mixed RF/FSO dual- 
hop transmission system. 



buildings. Since similar optical transmitters and detectors are 
used for FSO and fiber optics, similar bandwidth capabilities 
are achievable [18]. Therefore, this will get the required 
job done saving numerous amount of economic resources 
by utilizing FSO instead of digging up the current brown- 
field to install fiber optics between the different buildings. 
Another set of motivation behind such a system involves 
a fact that the users are mostly mobile and with only RF 
capabilities (no FSO capabilities). Installing FSO capability on 
these mobile users does not seem to be a justified approach. 
Simultaneously, we also fall short of bandwidth (BW) every 
now and then. Hence, to save on BW and to save on the 
economic resources by avoiding unnecessary modifications to 
the current mobile devices, we have introduced such a system 
wherein the users remain as is with RF only capability(s) 
and yet can be part of and/or make use of the FSO featured 
network. In another instance, we can think of a building floor 
(femto-cell in a heterogeneous network) where the users can 
send and receive through the backbone via FSO transmitter 
and detector respectively, placed at one of the corners of 
that floor. This FSO transmitter/detector can communicate 
with other such devices over other high-rise buildings and 
ultimately hop to the backbone. Additionally, to increase 
the spectral efficiency of such a system, we can study the 
effects on the performance of the system by selection of 
A -best users to be multiplexed. Hence, to perform such a 
study and due to space limitations, this manuscript tackles 
only the simplest possible scenario/special case of such a 
heterogeneous system and studies its statistical characteristics 
and ultimately the performance measures with on-going work 
on ultimately addressing all the issues mentioned above and/or 
as shown in Fig. 1 and beyond. 

However, the results presented in [1] were derived under the 
assumption of non-pointing errors in the FSO link but were 
limited to cumulative distribution function (CDF)/outage prob- 



We employ the same model as was employed in [1] and 
hence, the end-to-end SNR can be given as 7 = 71 7^ , 
where 71 represents the SNR of the RF hop i.e. S-R link, 
72 represents the SNR of the FSO hop i.e. R-D link, and C 
is a fixed relay gain [1], [10], [19]. 

The RF link (i.e. S-R link) is assumed to follow Rayleigh 
fading whose SNR follows an exponential distribution, pa- 
rameterized by the average SNR 7l of the S-R link, with 
a PDF given by / 71 (7i) = l/7i ex P(-7i/7i) l 19 l 0n 
the other hand, it is assumed that the FSO link (i.e. R- 
D link) experiences Gamma-Gamma fading with pointing 
error impairments whose SNR PDF is given under indirect 
modulation/direct detection (IM/DD) by [8, Eq. (12)], [9, Eq. 
(20)] that can be expressed in a simpler form by utilizing [20, 
Eq. (6.2.4)], as 



772(72) 



2 72 r( a )r(/3) 



ri3,0 
"1,3 



e 2 + i 
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where 7s is me average SNR of the R-D link, a and j3 are the 
fading parameters related to the atmospheric turbulence con- 
ditions [3]-[5] with lower values of a and /3 indicating severe 
atmospheric turbulence conditions, £ is the ratio between the 
equivalent beam radius at the receiver and the pointing error 
displacement standard deviation (jitter) at the receiver [8], [9], 
r(.) is the Gamma function as defined in [21, Eq. (8.310)], 
and G(.) is the Meijer's G function as defined in [21, Eq. 
(9.301)]. 



III. Closed-Form Statistical Characteristics 

A. Cumulative Distribution Function 
The CDF is given by [10] 



F 7 7 = Pr 



7i72 
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which can be written as 
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Using [22, Eq. (07.34.03.0228.01)], we can rewrite 
exp (-7 (7/(72 70) as Gj;?^-| n . Further using 



[20, Eq. (6.2.2)], 



W.O 



7 C 



72 7i 10 

, „ can be alternated to 

72 7i 1 J 

Now, along with the above modifications, 
applying [23, Eq. (21)] to (3), and some simple algebraic 
manipulations along with utilizing [20, Eq. (6.2.4)], the CDF 
of 7 can be shown to be given by 
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87r"r(a)r(/3) ' a " u ~ if^~ - For 1116 non-pointing errors case, 
when £ — > oo, it can be easily shown that the CDF in (4) 
converges to 



F 7 ( 7 ) = l-^ 2 cx P (-7/7 1 )G^ 



where k 3 = f , °d±, |, £±1, , and A 2 = 
agreement with [1, Eq. (15)]. 
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5. Probability Density Function 

Differentiating (4) with respect to 7 , using the product rule 
then utilizing [22, Eq. (07.34.20.0001.01)], we obtain after 
some algebraic manipulations the PDF in exact closed-form 
in terms of Meijer's G functions as 
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For the non-pointing errors case, when £ — » oo, it can be 
easily shown that the PDF in (6) converges to 
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C. Moment Generating Function 

The MGF defined as A4 7 (s) = E [e _7S ] can be expressed 
in terms of CDF as A4 7 (s) = s f Q e~ 7S F 1 (^)d r y. Using this 



equation by placing (4) into it and utilizing [21, Eq. (7.813.1)], 
we get after some manipulations the MGF of 7 as 
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When £ — > 00 (i.e. non-pointing error case), the MGF in (8) 
can be easily shown to converge to 
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D. Moments 



The moments defined as E [7™] can be expressed in terms 
of the complementary CDF (CCDF) Ffa) = 1 - F 7 ( 7 ) as 
E [7"] = n J °° 7 n_1 F^ (7)^7. Now, using this equation by 
placing (4) into it and utilizing [21, Eq. (7.813.1)], we get the 
moments as 



[7" 
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When £ — > 00, the moments in (10) can be easily shown to 
converge to 
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IV. Applications to the Performance of 
Asymmetric RF/FSO Relay Transmission Systems 

A. Higher- Order Amount of Fading 

The AF is an important measure for the performance of 
a wireless communication system as it can be utilized to 
parameterize the distribution of the SNR of the received signal. 
In particular, the n'^-order AF for the instantaneous SNR 7 
is defined as AF^ n) = E [ 7 ™]/E [7]" - 1 [24]. Now, utilizing 
this equation by substituting (10) into it, we get the n^-order 
AF as 

at?M A^-n n 6,i\ u 1 — n i K i 1 r,6,l \ u 0,/ci] 
AF) >=nA 1 G 2 ' 6 B G 2 ' 6 B -1. 

K 2 K 2 

(12) 

For n = 2, as a special case, we get the classical AF [25] as 

(13) 
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For the non-pointing errors case, when £ — > 00, it can be 
easily shown that the n t/l -order AF in (12) converges to 
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(14) 

For n — 2, as a special case, we get the classical AF [25], for 
non-pointing errors case, as 

AF = AF^ =2A 2 ~ 1 G\'\\b ^\g\'\\b °1 -1. 

' L K 3 J ' L K 3 

(15) 
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TABLE I 

BER Parameters of Binary Modulations 



Modulation 


P 


q 


Coherent Binary Frequency Shift Keying (CBFSK) 


0.5 


0.5 


Coherent Binary Phase Shift Keying (CBPSK) 


0.5 


1 


Non-Coherent Binary Frequency Shift Keying (NBFSK) 


1 


0.5 


Differential Binary Phase Shift Keying (DBPSK) 


1 


1 



B. Error Probability 

1) Average BER: Substituting (4) into [26, Eq. (12)] and 
utilizing [21, Eq. (7.813.1)], we get the average BER Pt of a 
variety of binary modulations as 
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where the parameters p and q account for different modu- 
lation schemes. For an extensive list of modulation schemes 
represented by these parameters, one may look into [26]-[29] 
or refer to Table I. For the non-pointing errors case, when 
£ — >• oo, the BER in (16) can be easily shown to converge to 
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2) Average SER: In [30], the conditional SER has been 
presented in a desirable form and utilized to obtain the average 
SER of M-AM, M-PSK, and M-QAM. For example, for M- 
PSK the average SER P s over generalized fading channels is 
given by [30, Eq. (41)]. Similarly, for M-AM and M-QAM, the 
average SER P s over generalized fading channels is given by 
[30, Eq. (45)] and [30, Eq. (48)] respectively. On substituting 
(8) into [30, Eq. (41)], [30, Eq. (45)], and [30, Eq. (48)], we 
can get the SER of M-PSK, M-AM, and M-QAM, as shown 
below 
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respectively. The analytical SER performance expressions ob- 
tained in (18), (19), and (20) are exact and can be easily esti- 
mated accurately by utilizing the Gauss-Chebyshev Quadrature 
(GCQ) formula [31, Eq. (25.4.39)] that converges rapidly, 
requiring only few terms for an accurate result [11]. 

C. Ergodic Capacity 

The ergodic channel capacity C defined as C = 
E[log 2 (l + 7)] can be expressed in terms of the CCDF of 
7 as C - l/ln(2)/ °°(l + 7 )- 1 / 1 7 c ( 7 )d 7 [32, Eq. (15)]. 
Utilizing this equation by exploiting the identity [33, p. 152] 



(l + az)- 
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b i,i 
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in it and using the integral 



identity [26, Eq. (20)], the ergodic capacity can be expressed in 
terms of the extended generalized bivariate Meijer's G function 
(EGBMGF) (see [26] and references therein) as 
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For the non-pointing errors case, when £ —> oo, the ergodic 
capacity in (21) can be easily shown to converge to 



— _ ^27l ^,1,0:1,1:5,0 



c 



ln(2) 



pi,u:i,i: o,u 
"1,0:1,1: 0,5 



1 



K3 



ll 



B 



(22) 



The expression in (21) and (22) can be easily and efficiently 
evaluated by utilizing the MATHEMATICA® implementation 
of the EGBMGF given in [26, Table II]. 

V. Results and Discussion 

The average BER performance of different digital binary 
modulation schemes are presented in Fig. 2 based on the values 
of p and q as presented in Table I. We can observe from 
Fig. 2 that the simulation results provide a perfect match to 
the analytical results obtained in this work. 

It can be seen from Fig. 2 that, as expected, CBPSK 
outperforms NBFSK. Also, the effect of pointing error can 
be observed in Fig. 2 i.e. as the effect of pointing error (as 
the value of £ increases, the effect of pointing error decreases) 
increases, the BER deteriorates and vice versa. It can be shown 
that as the atmospheric turbulence conditions get severe i.e. 
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Comparison between Analytical and Simulation Results 




Average Signal-to-Noise Ratio (SNR) per Hop (dB) 



VI. Concluding Remarks 

We derived novel exact closed-form expressions for the 
CDF, the PDF, the MGF, and the moments of an asymmetric 
dual-hop relay transmission system composed of both RF and 
FSO environments with pointing errors in terms of Meijer's G 
functions. Further, we derived analytical expressions for vari- 
ous performance metrics of an asymmetric dual-hop RF/FSO 
relay transmission system with pointing errors including the 
higher-order AF, error rate of a variety of modulation schemes, 
and the ergodic capacity in terms of Meijer's G functions. In 
addition, this work presents simulation examples to validate 
and illustrate the mathematical formulation developed in this 
work and to show the effect of the atmospheric turbulence and 
pointing error conditions severity and unbalance on the system 
performance. 
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